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FOREWORD 


This  Supplement,  prepared  under  Work  Unit  24041059,  Project  2404, 
in  the  Aerodynamics  and  Airframe  Branch,  Aerodynamics  Division,  Flight 
Dynamics  Laboratory  in  May  1982,  covers  improvements  in  the  method  of 
calculation  and  some  new  results  for  the  Technical  Report  AFWAL-TR-3105 , 
published  by  this  Laboratory  in  November  1981. 
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LIST  OF  SYMBOLS 


A  cross  sectional  area 

D  nozzle  diameter 

L  nozzle  segment  length  (also  nozzle 

length  in  the  calculation  results) 

M  flow  Mach  number 

p  static  pressure 

R  gas  constant 

t  segment  exit  to  inlet  area  ratio 

T  static  temperature 

v  gas  velocity 

cf  pipe  friction  coefficient 

Y  ratio  of  spec,  heats 

p  density 

T  fluid  force  parameter 

T  *  (t+i)/(2t),  see  Ref.  1 


INDICES: 

o  refers  to  stagnation  conditions 

1  refers  to  nozzle  segment  inlet 

2  refers  to  nozzle  segment  exit 

g  refers  to  gas 

p  refers  to  particle 
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I.  INTRODUCTION 


This  Supplement  to  Technical  Report  AFWAL-TR-8 1-3 105  (Ref.  1)  provides 
some  improvement  to  the  method  used  therein  for  calculating  the  flow  in  an 
expansion  nozzle  under  the  influence  of  wall  friction  and  heat  losses.  As 
mentioned  in  Reference  1  the  method  had  been  previously  developed  as  a 
general  tool  for  calculating  adiabatic  flows  in  ducts  with  friction.  By 
using  a  segmental  approach  the  method  was  also  considered  suitable  for 
calculating  non-adiabatic  flows  as  they  occurred  in  the  expansion  nozzle  of 
Reference  1.  In  this  case  the  nozzle  was  divided  into  many  small  segments. 
The  flow  in  an  individual  segment  was  then  treated  as  being  adiabatic  and 
only  between  segments  the  flow  conditions  were  corrected  for  the  occurring 
heat  loss.  This  process  has  now  been  found  to  be  a  poor  approximation  of 
the  real  flow.  New  investigations  revealed  that,  for  a  better  approxima¬ 
tion,  the  segments  themselves  must  be  treated  as  non-adiabatic ,  i.e.  the 
total  temperature  change  over  the  segment  must  be  part  of  the  input  to 
calculate  the  flow  in  the  segment.  Since  the  total  temperature  change  is  a 
function  of  the  unknown  flow  conditions  in  the  segment,  an  iteration 
procedure  appears  to  be  necessary.  However,  calculations  have  shown  that 
for  small  enough  steps  for  which  the  heat  transfer  conditions  do  not  change 
abruptly  from  segment  to  segment,  it  is  sufficient  to  simply  assign  the 
total  temperature  change  of  the  preceeding  segment  to  the  following  one  for 
calculating  the  non-adiabatic  flow  conditions.  Heat  loss  and  total  temper¬ 
ature  change  in  the  segment  are  then  recalculated  before  proceeding  to  the 
next  segment . 

The  improved  method  has  been  checked  for  accuracy  by  comparing  its 
results  with  those  obtained  for  the  same  case  with  Shapiro's  method  of 
calculating  flows  in  ducts  described  in  Reference  2.  For  each  method  the 
steps  in  the  calculation  procedure  were  chosen  small  enough  to  practically 
eliminate  the  influence  of  the  step  size  on  the  results.  Under  these 
conditions  nearly  perfect  agreement  between  the  two  methods  was  obtained. 
The  comparison  of  these  two  methods  together  with  discussions  on  the 
practicability  of  each  method  has  been  made  a  subject  of  a  separate  report 
intended  for  publication  as  an  AIAA  Technical  Note  (Ref.  3). 
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The  following  sections  contain  the  derivations  for  the  improved 
calculation  method  and  some  new  results  for  the  various  nozzle  cases  of 
Reference  1. 


II.  THE  IMPROVED  CALCULATION  METHOD 

Only  Equation  14  in  Reference  1  is  affected  by  the  new  calculation 
method.  To  explain  the  change,  the  derivations  given  in  Reference  1  are 
partly  repeated  here  beginning  with  Equation  10,  with  the  changes  intro¬ 
duced  which  lead  to  an  improved  Equation  14.  To  better  understand  the 
changes,  the  derivations  have  been  made  more  detailed. 


+ 


1) 


(1) 


This  equation,  derived  from  the  flow  momentum  conditions,  contains  two 
unknowns:  the  exit  Mach  number  M^  and  the  pressure  ratio  p2  /p  j  across 
the  expanding  flow.  Conservation  of  mass  and  energy  as  expressed  by  the 
following  two  relations 

V1  *1  C1  "  V2  k2  °2  (2)  . 


(3) 


allows  one  to  write  an  additional  condition  for  the  two  unknowns.  The 
energy  equation  under  the  assumption  of  ideal  gas  conditions*) ,  describes 

*)  The  ideal  gas  assumption  is  very  realistic  for  helium  used  in  the 
present  calculations.  For  gases  with  state-dependent  gas  properties  too, 
the  ideal  gas  assumption  constitutes  no  severe  restriction.  In  the 
segmental  approach  the  assumption  applies  only  to  a  small  segment  of  the 
expansion  nozzle  and  the  gas  properties  can  be  changed  from  segment  to 
segment.  For  a  rougher  approximation  constant  average  gas  properties  may 
bo  assumed  for  the  whole  expansion  process  as  done  for  the  case  of  air  in 
Table  1. 
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lated  to  the  total  temperature  by  the  energy  equation  (Equation  3) . 


Let 


t 


then  Equation  3,  9,  and  10  combined  give 


m2(1  +  ^  m22>1/2 


^  Pi  l/<V 

c  p2  Y  Ti 


(Vo 

'Vo 


(ID 


This  equation  is  the  "improved  equation  14".  The  change  from  the 
original  equation  consists  of  the  addition  of  the  ratio  of  the  segment 
exit-to-inlet  total  temperature.  This  ratio  must  be  derived  from  the  heat 
transfer  occurring  on  the  segment  walls.  This  heat  transfer  is  a  function 
of  the  still  to  be  determined  flow  conditions  in  the  segment.  To  avoid  an 
iteration  process  in  calculating  heat  transfer  and  flow  conditions,  it  is 
sufficient  to  use  the  total  temperature  ratio  of  the  preceeding  segment  to 
calculate  the  flow  conditions  in  the  new  segment  by  means  of  Eq .  11.  For 
the  first  segment  it  is  sufficient  to  use  unity  for  the  total  temperature 
ratio,  i.e.  adiabatic  conditions  as  done  in  Ref.  1. 


III.  RESULTS 

The  subsonic  nozzle  flow  remains  practically  unaffected  by  the  im¬ 
proved  method  due  to  the  small  change  in  the  total  gas  temperature  occur¬ 
ring  in  this  region.  For  easy  reference,  the  old  Figure  6  showing  the 
subsonic  conditions  is  included  here  as  Figure  1.  The  notations  on  this 
figure  also  provide  the  unchanged  input  data  for  the  calculations.  The 
data  list  has  been  extended  to  indicate  that  the  pipe  friction  coefficient 
is  0.02  for  the  subsonic  and  0.015  for  the  supersonic  nozzle  portion. 
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Figures  7a  -  c  of  AFWAL-TR-81-3105  have  been  recalculated  with  the  improved 
method.  The  results  are  shown  here  as  Figures  2a  -  c.  Some  of  the  less 
essential  curves  have  been  omitted  from  the  figures  to  improve  their 
clarity. 

To  simplify  the  calculations,  the  new  Figures  2b  and  2c  are  calculated 
without  the  short  0.34°  nozzle  portion  next  to  the  throat,  which  was 
included  in  the  old  calculations  (Figures  7b  and  7c).  This  0.34°  portion 
is  the  result  of  a  design  consideration  and  not  essential  for  the  particle 
acceleration.  The  nozzle  portion  near  the  throat  requires  a  special  design 
for  handling  the  enormous  heat  load  in  this  region.  For  actual  experi¬ 
ments,  extensions  with  various  wall  angles  would  be  attached  to  this 
specially  designed  throat  portion. 

The  improved  method  affects  mostly  the  supersonic  Mach  numbers, 
increasing,  for  instance,  in  the  case  of  the  0.34°  nozzle,  the  maximum  Mach 
number  from  1.4  to  2.05.  It  also  strongly  increased  the  drop  in  the  static 
gas  pressure  along  the  nozzle.  Since  the  maximum  gas  velocity  was  only 
increased  from  3280  m/sec  to  3730  m/sec,  the  large  increase  in  Mach  number 
must  be  associated  with  a  substantial  decrease  in  the  local  velocity  of 
sound  or  under  the  present  circumstances,  with  a  substantial  decrease  in 
the  static  gas  temperatures.  Since  the  heat  losses  were  reduced,  i.e.  the 
drop  in  total  gas  temperature  was  decreased,  the  lowering  of  the  static 
temperatures  is  essentially  the  consequence  of  the  faster  static  pressure 
drop  indicated  above.  Since,  as  also  indicated,  the  gas  velocities  were 
not  greatly  changed,  the  particle  velocities  were  even  less  affected, 
changing  in  the  present  example  from  1950  m/sec  to  2100  m/sec  for  the 
maximum  value  obtained. 

In  the  original  calculations,  a  nozzle  wall  angle  of  0.5°  was  found  to 

be  optimum  for  the  acceleration  of  the  particles.  The  new  calculations 

shift  the  optimum  to  smaller  angles.  As  an  example,  the  performance  for  a 

0.25°  wall  angle  is  shown  in  the  new  Figure  2d.  The  performance  is  also 

listed  in  Table  I  as  No.  3a.  For  p/p  *  0.02,  this  nozzle  gives  the 

o 

highest  particle  velocity  for  all  nozzles  with  a  throat  diameter  of  3.97mm. 
'.he  diffeT  nee  in  particle  velocities  between  these  four  nozzles  is  not 
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very  pronounced.  Considerations  for  the  design  and  fabrication  of  nozzles 
with  such  small  wall  angles  may  dictate  the  actual  choice.  Also,  the 
nozzle  length  may  be  a  decisive  factor.  With  a  small  sacrifice  in  particle 
velocity,  fairly  short  nozzles  with  a  large  wall  angle  give  good  results. 
For  instance,  case  No.  2  in  Table  1  assumes  a  wall  angle  of  0.5°  and 
results  in  a  nozzle  length  of  only  344  mm  to  achieve  a  particle  velocity  of 
2014  m/sec. 

A  note  for  Table  II  of  AFWAL-TR-81-3105:  the  indication  of  a  minimum 
wall  angle  in  case  No.  9  was  the  result  of  misinterpreted  computer  results. 
The  given  value  is  incorrect.  No  attempt  has  been  made  to  find  the  correct 
value.  In  concurrence  with  the  above  discussion  on  optimum  wall  angles, 
one  can  only  conclude  that  for  the  present  wall  friction  assumptions  the 
minimum  wall  angles  for  which  the  flow  becomes  choked  are  in  all  cases 
extremely  small  with  case  No.  9  being  no  exception. 
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CALCULATION  RESULTS 
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Figure  1.  Flow  Conditions,  Heat  Loss,  and  Particle  Acceleration  for  the 
Subsonic  Portion  of  the  Expansion  Nozzles  Considered  in 
Figures  2a  to  d  for  Their  Supersonic  Nozzle  Portions 


Flow  Mach  number  M,  Gas  Velocity  V-Gas,  Particle  Velocity  V-Part 
Total  Gas  Temperature  T0,  Expansion  Pressure  Ratio  P/PG  and 
Accumulated  Heat  Loss  from  the  Nozzle  Q-Loss  for  the  Supersonic 
Nozzle  Portion  at  Different  Wall  Angles  (friction  coefficient  Cf 
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